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Multi-Antenna Wireless Energy Transfer for 
Backscatter Communication Systems 
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Abstract 

We study RF-enabled wireless energy transfer (WET) via energy beamforming, from a multi-antenna 
energy transmitter (ET) to multiple energy receivers (ERs) in a backscatter communication system, such 
as REID, where each ER (or REID tag) reflects back a portion of the incident signal to the ET (or REID 
reader). Eor such a system, the acquisition of the forward-channel (i.e., ET-to-ER) state information (E-CSl) 
at the ET is challenging, since the ERs are typically too energy-and-hardware-constrained to estimate or 
feed back the E-CSl. The ET leverages its observed backscatter signals to estimate the backscatter-channel 
(i.e., ET-to-ER-to-ET) state information (BS-CSI) directly. We first analyze the harvested energy obtained by 
using the estimated BS-CSI. Furthermore, we optimize the channel-training energy and the energy allocation 
weights for different energy beams, for weighted-sum-energy (WSE) maximization and proportional-fair- 
energy (PEE) maximization. For WET to single ER, we obtain the optimal channel-training energy in a 
semi-closed form. For WET to multiple ERs, the optimal WET scheme for WSE maximization is shown 
to use only one energy beam. For PEE maximization, we show it is a biconvex problem, and propose a 
block-coordinate-descent based algorithm to And the close-to-optimal solution. Numerical results show that 
with the optimized solutions, the harvested energy suffers slight reduction of less than 10%, compared to 
that obtained by using the perfect F-CSI. Hence, energy beamforming by using the estimated BS-CSI is 
promising, as the complexity and energy requirement is shifted from the ERs to the ET. 

Index Terms 

Backscatter communication systems, wireless energy transfer, energy beamforming, channel estimation, 
resource allocation, proportional fairness, biconvex optimization 

I. Introduction 

Recently, backscatter radio has been utilized widely, due to its low energy requirement and low 
monetary cost. Backscatter radio performs communication by means of reflection of incident signals 
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rather than direct radiation, The most prominent commercial use of backscatter radio is in radio 
frequency identification (RFID) applications, for identifying people or products in supply chains. For 
backscatter communication, the input impedance of the tag’s antenna is intentionally mismatched to 
scatter back a portion of the incident signal. The phase and amplitude of the backscattered signal is 
then determined by the input impedance. By varying the antenna impedance, an RFID tag encodes 
digital symbols into the backscattered signal, which is then received and decoded by the RFID 
reader m. The tag operates without any on-tag energy source and relies entirely on backscatter, 
which leads to its energy efficiency and cost effectiveness. The existing literature on communication 
theoretic aspect of backscatter focuses on the tag-to-reader channels, such as on multiple access 
techniques [|3, performance of space-time code as well as the diversity-multiplexing tradeoff 
for multi-input-multi-output (MIMO) [|T1|. 

Backscatter communication is also valuable to other systems built with low-power and low-cost 
principles, such as wireless sensor networks. The ongoing integration of various sensors on RFID 
tags, such as described in ||5l lEl, confirms the potential for new sensor networks that use modified 
RFID components to transfer sensor data to the fusion center. An example of an RFID-based sensor 
platform is the wireless integrated sensing platform (WISP) Q. 

However, the coverage range of RFID sensors is significantly limited by the forward channels 
(i.e., reader-to-tag), as only a small amount of RF energy can be harvested at the conventional 
RFID sensors. In the literature, the range of commercial RFID tags is improved from aspects 
including rectifier circuit design (see 0 and references therein) and special waveform design [[8]|. 
The waveform is designed such that the received signal frequently exceeds the threshold voltage 
required to turn on the rectifier circuit. 

Multi-antenna techniques have been shown to be efficient for enhancing the efficiency of wireless 
energy transfer (WET) for traditional radio communication systems. The electromagnetic (EM) 
energy needs to be concentrated into a narrow beam to achieve efficient transmission of energy, 
referred to as energy beamforming dH, as EM waves decay quickly over distances. The channel 
state information (CSI) is prerequisite for energy beamforming. The ERs perform channel estimation 
(CE) by receiving pilots sent from the FT, and then feed back the estimated forwarcii] channel (i.e., 
ET-to-ER) state information (F-CSI) to the FT; or the FT receives pilots sent from the ERs, and 


*For consistence, the terminology of “downlink channel” and “uplink channel” in traditional radio communication systems is 
renamed to he “forward channel” and “backward channel”, respectively. 
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then obtains the estimated F-CSI directly by exploiting the channel reciprocity. The effect of CE 
and feedback on energy beamforming was studied in [fT^ - lfT3ll . In particular, ifTOTl investigated the 
dynamic allocation of time resource for CE and energy resource for WET. ffTTl studied energy 
beamforming by using one-bit feedback from each ER, to facilitate hardware implementation. 
Based on the one-bit feedback, the ET adjusts transmit beamforming and concurrently obtains 
improved estimates of the forward channels to all ERs. ifT^ maximized the net harvested energy 
after subtracting the energy used for the ER sending pilots, for a point-to-point MIMO WET system. 
Eurthermore, with the estimated E-CSI, lfT3l optimized the throughput for a massive MIMO system 
powered by WET. 

Eor backscatter communication systems, energy beamforming can also be used to achieve efficient 
WET. Multiple antennas are deployed at the ET (or REID reader) to perform energy beamforming 
toward ERs (or REID tags). However, unlike traditional radio communication systems, the acquisition 
of the E-CSI at the ET is challenging, since the ERs are typically too energy-constrained to perform 
CE or feedback, and also may not have specific hardware built for CE nor feedback. 

Instead, the ET may in fact leverage its observed backscatter signal to estimate the backscatter- 
channel (i.e., ET-to-ER-to-ET) state information (BS-CSI) solely by itself. This shifts complexity 
and energy requirements from the ERs to the ET. The necessary synchronization condition for CE 
can also be achieved more easily for the ET than for the ERs. A proof of concept was presented 
in m, which illustrated that the WET can be optimized by using only the power levels received 
at the receive antennas of the ET. However, in llT4l . the ET transfers energy to only one ER in 
a time period, even though other ERs may desire energy and can also potentially harvest energy 
concurrently. Moreover, that work did not consider the effect of CE nor the resource allocation for 
WET to multiple ERs concurrently. 

In this paper, based on a backscatter communication system, we consider the WET from an 
ET with multiple antennas to multiple ERs each with a single antenna. We assume frame-based 
transmissions, where each frame consists of a CE phase and a WET phase. With unknown backward 
channels (i.e., ER-to-ET), we first analyze the energy that can be harvested via energy beamforming 
by using the estimated BS-CSI. Eurthermore, we optimize the resource allocation, by investigating 
two utility maximization problems, namely, the weighted-sum-energy (WSE) maximization and the 
proportional-fair-energy (PEE) maximization. The optimization variables are the channel-training 
energy used for CE, and the energy allocation weights for beamforming toward multiple ERs. Eor 
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comparison, we consider two benchmarks, one is the ideal case of energy beamforming by using the 
perfect F-CSI, in which the most energy can be transferred via energy beamforming in any wireless 
communication system; the other is energy beamforming by using the estimated F-CSI fed back 
from the ERs, in a traditional radio communication system. 

The main contributions of this paper are as follows: 

• We propose a novel energy beamforming scheme by using the estimated BS-CSI, to perform 
WET to multiple ERs concurrently. This scheme shifts the complexity and energy requirements 
from the ERs to the ET, and is thus especially attractive for transferring energy to (ultra-)low- 
power and low-cost wireless devices that can neither estimate channels nor send pilots or 
feedback. 

• We obtain an analytical expression for the harvested energy obtained by using the estimated 
BS-CSI, in which the ambiguity of unknown backward channels is taken into account. We also 
obtain bounds on the harvested energy, which are numerically shown to be tight. 

• We obtain the optimal resource allocation schemes for both WSE maximization and PEE 
maximization. Eor the single-ER case, we obtain the optimal channel-training energy in a 
semi-closed form. Eor the multiple-ER case, to achieve WSE maximization, the optimal WET 
scheme is shown to use only one energy beam. Eor PEE maximization, we show it is a biconvex 
problem, and propose a block-coordinate-descent (BCD) based algorithm to find the close-to- 
optimal solution. Numerical results show that fairness is improved by PEE maximization. 

• We conduct simulation studies, which show that the maximally harvested energy suffers slight 
reduction of less than 10% when compared to the harvested energy obtained by using the 
perfect E-CSI, and of about 3% when compared to the net harvested energy obtained by using 
the estimated E-CSI in a traditional radio communication system. Eor the latter benchmark, 
the net harvested energy is the harvested energy after subtracting the energy used for sending 
backward pilots and feeding back the estimated E-CSI. This observation is encouraging, as 
almost all complexity of the ER is shifted to the ET. 

The rest of this paper is organized as follows: Section UH presents the system model. Section 
rrm analyzes the harvested energy. Section UV] considers the resource allocation, by formulating a 
general utility-maximization problem. Section |V] obtains the optimal resource allocation for WSE 
maximization and PEE maximization. Section gives some discussion. Section IVIII provides 
extensive numerical results. Einally, Section IVIIII concludes this paper. 
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II. System Model 

As illustrated in Fig. [H based on a backscatter communication system, we consider WET from an 
ET (or RFID reader) that can concurrently transmit with M antennas and receive with R antennas, 
to K ERs (or RFID tags) each with single antenna. The ET can perform CE, energy beamforming 
and other signal processing operations. Each ER contains an RE-energy harvesting module which 
supplies energy for operations such as sensing, quantization and backscatter modulation. Besides, it 
has a switched load impedance that is connected to its antenna. By varying the antenna impedance, 
each ER encodes digital symbols into the backscattered signal, which is then received and decoded 
by the ER. Compared to traditional radio communication systems, backscatter communication is 
energy efficient, as it just reflects the incident signal, without generating radio signals actively. 
Hence, backscatter communication is alluring for (ultra-)low-power devices Q. 



Eig. 1: System model 


We study WET via frame-based transmissions on a single frequency band. As illustrated in 
Eig. [2l the time duration of each frame is fixed as T symbol periods, which consists of the CE 
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phase followed by the WET phase. The CE phase of r symbol periods is further equally divided 
into K slots, eaeh of whieh eonsists of ML sueeessive symbol periods. That is, the CE time is 
r = KML. During the CE phase, the ET sends pilot signal with power pi. With the eoordinatiorO 
of the ET, in the k-th CE slot, only the load impedance of the k-th ER is switched on to facilitate 
backscatter communications, with the impedance of all other ERs switched off. The pilot signal is 
thus backscattered by only the k-th ER. After receiving the backscatter signal, the ET estimates the 
backscatter-channel associated to the k-th ER. During the WET phase of (T — t) symbol periods, the 
ET performs energy beamforming, and all ERs switch off the load impedance and harvest wireless 
energy. Typically, the time duration of the CE phase is much shorter than the WET phase. Hence, 
we assume that the ERs do not harvest energy during the CE phase, for simplicity. 

A. Backscatter Channel 

The backscatter channel is modeled as a concatenation of three components, namely, forward 
channel (i.e., reader-to-tag), backscatter reflection coefficient, and backward channel (i.e., tag-to- 
reader). Eet hmk denote the forward channel between the m-th transmit antenna and the k-th ER, 
and gkr denote the backward channel between the k-th ER and the r-th receive antenna. Denote the 
(long-term) path loss of the channel between the ET and k-th ER by (3k, which is assumed to be 
constant over frames and taken to be known a priori at the ET. We assume the forward and backward 
channels are flat Rayleigh-fading and independent, i.e., the channel coefficients hmk ~ CA/’(0,/9fc) 
and Qkr ~ CAA(0,/3fc). In the k-th CE slot, the load impedance of ER k is switched on, and ER 
k reflects a portion of the incident signal to the ET, which is modeled by the complex reflection 
coefficient pk G C. Eor other ER j, the impedance is switched off, i.e., pj = 0, Vj ^ k. In this 
paper, we assume the ERs are linear backscatter devices, i.e., the reflection coefficient pk is fixed 
and does not vary with the incident power at the ERs. Without loss of generality, we assume pk = 1. 

The backscatter channel associated with the m-th transmit antenna, the k-th ER and the r-th 
receive antenna is given by Umkr = hmkQkr- Hence, the backscatter channel experiences double 
fading due to hmk and pkr- Eor convenience of expression, we let the number of receive antenna 
R = 1 and omit the subscript r in the notation. The analysis can be easily extended to the case of 
multiple receive antennas, but is beyond the scope of this paper. 

^The ET can coordinate all ERs to switch on in any arbitrary sequence. For instance, the ET sends an initialization symbol to 
signal the start of the CE phase, and then transmits a unique ID at the beginning of each CE slot. The corresponding ER then 
responds accordingly. 
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Frame 1 


Channel Estimation 


Frame 2 


Frame ... 


Wireless Energy Transfer 
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Fig. 2: Frame structure 


B. Channel Estimation via Backscatter Signal 


The total energy for sending pilots during the whole CE phase is fixed as rpi. We assume the ET 
spend equal pilot energy (i.e., MLpi) for estimating the backscatter channel associated with each 
ERy. Under the assumption of independent channels, from tfTSll . the least-square (ES) estimation 
performance can be optimized by using the pilot matrix X = [Xi X 2 ■ ■ ■ X/,]^ in each slot, where 
X; is an orthogonal matrix such as the Hadamard matrix or the identity matrix, with power pi, 
i.e., X/Xf^ = piIm- The (■)^ and (■)^ denote transpose and conjugate transpose, respectively. Eor 
convenience, we choose X/ = ^/p^M, I = - ,L. 

When the ET sends pilots, the pilot signal can also be received directly by its own receive antenna, 
without going through the backscatter channel. Eor simplicity, we assume this undesired signal can 
be estimated accurately and subtracted from the received signal, as this direct channel is static and 
can be estimated a priori. Hence, the received signal with respect to the pilot sent from the m-th 
transmit antenna is given by 

L 

Vink T 'y ^ '^mkl 1 (1) 

l=l 


where the ET noise rimki^ are independent and distributed as CAA(0,(J^). The ES estimate for the 
backscatter channel with respect to the m-th transmit antenna and the fc-th ER is obtained as 


a^k = = a^k + 


L^/^ 


( 2 ) 


^In general, the training energy for different ERs can be adjusted depending, for example, on the long-term path loss. However, 
that is beyond the scope of this paper. 
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where the equivalent noise ~ CN'{Q, -f^). For energy beamforming, the estimate of the F-CSI 


w ~ ^u, 

hrak is desired. With the estimated BS-CSI amk, there is however still a remaining ambiguity due to 


the unknown baekward ehannel gk- We will see later the effeet of the ambiguity Qk on the harvested 
energy, in Seetion ITTT-AI and Seetion ITTT-Bi 

However, the performanee of CE depends strongly on g^. We define the intermediate random 


variable hmk — ^mkgk the error e-mk — hmk — hmk- Conditioned on gk, it is then standard to 


show that the forward ehannel hmk is eonditionally distributed as 

l^khmk 


hmk\(^mk, gk ~ CJ\f 


fk + (^lki9k)'Pk + (^lki9k)^ 

and the error is eonditionally distributed as Cmk ~ CAA(0, (Tl^{gk)), with error varianee 

KMa^ 


(3) 


^e,k{9k) = 


\9k\‘^Tpi' 


(4) 


Remark 1 (Effeet of unknown backward channel on forward-channel estimation). Eor the special case 
of = 1 and hence amk = hmk, the estimated E-CSI hmk is given in (O, which is exactly the same 
as that for a traditional radio communication system. However, for a backscatter communication 
system, from dH), the estimation error depends on the (unknown) backward channel gk, as the signal- 
to-noise (SNR) for CE is affected by gk. Thus we have to account for the unknown g^ to derive the 
statistics of the expected harvested energy later in Section |nll 


HI. Wireless Energy Transfer via Energy Beamforming 

In this section, we study the WET by using the estimated BS-CSI. In order to analyze how it 
differs from WET by using the estimated E-CSI in traditional radio communication systems, we first 
consider the single-ER case in Section IIII-Al and analyze the effect of unknown backward channel 
on the harvested energy in Section IIII-Bl The case of multiple ER is studied in Section Illl-Ci in 
which we obtain bounds on the harvested energy, so as to simplify analysis in next sections. 


A. WET to Single ER 

Instead of considering the specific case of WET to one ER, for full generality, we consider the 
case of WET to ER k. Eor beamforming toward only ER k, the M x 1 transmitted signal is given 
by y'^w(afc), where p 2 is the transmit power for WET, and w(afc) is the beamformer depending 
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on the estimated BS-CSI a^. We note that traditionally the estimated F-CSI is used to obtain the 
beamformer. The reeeived signal by ER k is written as 

Zk = vSw^(afc)hfc + Hk. (5) 


where the noise at ER k is distributed as CAA(0,cr^Q). Due to the law of energy eonservation with 
effieieney rj, the RF-band energy harvested by ER k during the WET phase, denoted by Ek, is 
assumed to be proportional to that of the reeeived baseband signal, i.e.. 


Ek = Ek{T,pi,p2) = pP2{T - r)Eh^ 




W 


(^k)hk 


( 6 ) 


We assume in db]) that the energy due to the noise at ER k eannot be harvested. Eor eonvenienee, 
we also assume p = I'm this paper. Eet (•)* denote the eomplex eonjugate. The expeeted harvested 
energy in ® is rewritten as 


Ek{,T,Pi,P2) =P2 {T-t)E, 


9k 


E- 


\9k 


w 


(a) 


P 2 {T - r)E 


'9k 


E- 


'SfelSfe 


afc)Eh,|s„g, [hfch^] w*(afc) 

l3k(Tlki9k)lM , /Sfchfchf 


w'^(afc 


+ 


+ (Ze,ki9k) [/3k + (Tlk{gk)y 


w*(afc) 


, (V) 


where the eonditional eorrelation matrix in (a) is obtained from the eonditional distribution in 
dS]). Unlike traditional radio eommunieation systems, the harvested energy in dU) is obtained by 
performing expeetation over Pk, As noted in Remark [H this takes into aeeount the ambiguity of pk- 
Net, we obtain the optimal beamformer (with unknown pk) in the following lemma. 


Lemma 1. The optimal energy beamformer toward a single ER k is given by 

w(afc) = 

11 11 

The corresponding maximally harvested energy by ER k is given by 


( 8 ) 


Ek{r,Pi,P 2 ) = P2M(3k{T - r) I 1 - , ^ E 


M 


'9k 


3kTpi\9k\^ 

L KMcr^ 


+ 1 


(9) 


Proof: See Appendix 


Remark 2 (Diseussion for the speeial ease of large number of transmit antennas M or users K). As 
the produet MK tends to infinity, the harvested energy in dH) approaehes the quantity P 2 Pk{T — r). 
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This is equivalent to the ease when the wireless energy is harvested from the omnidireetional signal 
transmitted by the ET. This ean be explained intuitively as follows. For large M or K, a finite 
amount of energy has to be shared for training over all ehannels between the transmit antennas and 
the antennas of all ERs, whieh leads to inaeeurate estimates of the BS-CSI. Henee, the beamforming 
gain that ean be aehieved for WET is very limited. 

B. Effect of Unknown Backward Channel on Harvested Energy 

In this seetion, we analyze the effeet of unknown baekward ehannel gk on the harvested energy, 
still assuming WET to single ER for exposure. For radiative eommunieation systems, the optimal 
beamformer is the normalized estimated F-CSI [fTOl . From EemmalU for baekseatter eommunieation 
systems in whieh the baekward ehannel gk is unknown, the optimal beamformer that aehieves the 
maximally harvested energy in dH) is the normalized estimated BS-CSI. In the proof for Eemma [U 
we show that when gk is given, the optimal beamformer is just the normalized estimated F-CSI. 
Intuitively, this is beeause gk is eommon for estimating all the forward ehannels between the ET 
and the ER. Similar observation was also obtained in [|T4l . 

However, this ambiguity of baekward ehannel gk results in a reduetion of harvested energy, as 
shown the following proposition. 

Proposition 1. Assuming the same reeeive SNR for CE, the harvested energy Efc(r,pi,p 2 ) for a 
baekseatter eommunieation system is upper bounded by that for a traditional radio eommunieation 
system, given by 



( 10 ) 


Proof: Define the random variable Yk = Is'fcp, whieh follows exponential distribution. We 


further define the funetion f{yk) = cky\+i ^ where Ck = jdkrpi/{KMa"^). It ean be easily eheeked 
that f{yk) is a strietly eonvex funetion of yk for yk > 0, henee, EyJ/(lfc)] > 

^Yk [^fc] = Pk, due to Jensen’s inequality. Thus, the harvested energy in ® is upper bounded as in 
(QOl). This proves Efc(r,pi,p 2 ) < Efc(r,pi,p 2 )- 

For a traditional radio eommunieation system, when the reeeive SNR for CE at the ER is fdkCk 
whieh is the same as that for a baekseatter eommunieation system, it ean be shown that the harvested 
energy is exaetly Ek{T,pi,p 2 ) in (fTO . by following the steps in ifTOl . This eompletes the proof. ■ 
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Numerical results in Section IVIII will show that this reduction in the harvested energy is marginal, 
which motivates the use of backscatter WET due to its low complexity at the ER. 


C. WET to Multiple ERs 

To achieve WET to all ERs concurrently, we allow the use of multiple energy beams each 
toward one particular ER. Denote A = [ai,a 2 , ••• ,ai^]. The beamformer is then chosen as a 
linear combination of the normalized estimated BS-CSI a^’s, i.e., 

w(A) = W (11) 

where the weights ^^’s are subject to the condition Yl!k=i = 1- 

Remark 3. As the number of antennas M tends to infinity, the beamformer in (fTTl) is asymptot¬ 
ically optimal [fT3]l . This is the important motivation for choosing the beamformer in (fTTl) . while 
maintaining the flexibility of choosing different weights for energy beams toward different ERs. 


Similar to the case of single ER, the energy harvested by ER k is given by 


Ek = Ek{T,pi,p2,^k) = P2{T - r)E 


hfc,A 


w'^(A)h^ 


( 12 ) 


The harvested energy is given by the following lemma. 

Lemma 2. With the beamformer in (fTTl) . the harvested energy by ER k is given by 

Ek{,T,pi,p2,^k) = P2M(3kik{T- t) 

Proof: See Appendix fBf ■ 

The first term in (fT3l) is the harvested energy from the beam directly toward ER k, while the 
second term is the energy harvested from beams toward other ERs but still harvested by ER k. 

Erom Eemma 2 in lfT3ll . the harvested energy of a traditional radio communication system is the 
same as in (fT3l) for = 1- Hence, the backscatter communication system studied here generalizes 
the result of a traditional radio communication system. 

We note that the harvested energy in (fT3l) appears analytically intractable. Therefore, we will 
obtain bounds for the harvested energy, to simplify analysis in subsequent sections. Before that, we 
give the following lemma. 


A-rpilgfcP I ^ 


. KMa'^ 


+P2/3fc(T-r)(l-^fc). (13) 
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Lemma 3. Let random variable X ~ CAA(0,/3) and c be some positive eonstant. Then, 


Ex 


1 

1 + C|X|2 



(14) 


where Ti{t) = u ^exp(—M)dM is an upper ineomplete Gamma funetion. Moreover, the expee- 
tation is lower and upper bounded as 


ln(l + 2c/3) 


< E 


X 


l + c|X| 


< 


ln(l + c/3) 

c (3 


(15) 


Proof: The expeetation in (fT4l) is obtained by standard integration. The lower bound and upper 
bound in (fTSl) is obtained from [fT^ (5.1.20)]. ■ 

We assume the average power for eaeh frame is Pave- Besides the energy eonsumption for ehannel 
training, all the remaining energy is used for WET. This implies the WET power is given as 


P2 = 


PaveT - rpi 

T-t 


(16) 


Erom (fT4l) and (fT^ . the harvested energy in (l9|) is rewritten as 


Ek{r,Pi,^k) = Mf3k^k{PayeT - rpi) 


f3k{PaveT - rpi)(l - ^fc) 


M-lKMa"^ fKMa^\ 


M 


f^lrpi V Plrpi J ^ V /^irp, ) 


+ 


= PkiPaveT - q) 
— Ek{q,ik), 


(M-1) 


KMa"^ fKMa‘^\^ f KMa'^ 


V Ph / \ ah 


■Cfc + 1 


(17) 


where in (a) we perform algebraie simplifieation and introduee q = rpi whieh denotes the total 
energy used for ehannel training. We observe that the harvested energy depends only on q, sinee the 
CE time r and the transmit power pi for CE are always eoupled in (fTTI) . Based on this observation, 
we will optimize the ehannel-training energy q in next seetions, instead of optimizing r and pi. 

Eor the speeial ease that g = 0, no CSI is available and henee the ET performs omnidireetional 
transmission; then the harvested energy reduees to Ek{0,^k) = PaveT/Sk, as expeeted. 

Erom Eemma[3]and (fTTI) . we immediately obtain the bounds on the harvested energy as follows: 
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Lemma 4. With the energy beamformer in (fTTl) . the energy harvested by ER k is lower bounded by 

KMa\ A Pig V 


Ek{.q,ik) > f^kiPaveT - q) 


(M- 1) 


1 - 


„ In ( 1 “t" , 

I3lq V KMa^J 


Cfc + 1 


^Efc(g,eA:), (18) 


and it is upper bounded by 


Ek{q,^k) < l3k{p aveT 


q) 


(M-l) 



KMa^ 


In 



Wlq Y 

KMa^ ) _ 


+ 1 


(19) 


As will be numerieally shown in Seetion IVIIi the bounds in (flSl) and (fT9l) are tight, espeeially 
for the lower bound. In the sequel, for analytieal traetability, we take the lower bound Ek{q,^k) as 
the energy harvested by ER k. 


IV. Resource Allocation for WET in Backscatter Communication Systems 

In this seetion, we further eonsider the resouree alloeation for WET in a baekseatter eommunioa- 
tion system. We formulate a general optimization problem to maximize a total utility of the energy 
harvested by all ERs. Eet fkiEk) be the utility of the k-th ER’s harvested energy Ek given by (fT^ . 
assumed to be a monotonieally inereasing funetion of Ek. Denote the veetor of energy alloeation 
weights for different energy beams by ^ = [^i ,^2 • • • The total utility is then given by 

K 

U{q,^) = Y,h{Ek{q,^k)). (20) 

k=l 

We aim to maximize the total utility by optimizing the ehannel-training energy q and the energy 
alloeation weights ^ for different energy beams, subjeet to the total energy eonstraint (i.e., the training 
energy q ean not exeeed the total energy PaveT available in eaeh frame) and the normalization 
eonstraint for the energy alloeation weights The general utility maximization problem is thus 
formulated as follows 


K 

(Pi) max '^fk (Ek{q,^k)^ 

^ k=i 

(21a) 

S. t. 0 < g < PaveT 

(21b) 

K 



(21e) 

k=l 


6 > 0, Wk. 

(21d) 
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In particular, we will investigate two utility maximization problems, namely, the weighted-sum- 
energy (WSE) maximization where the utility funetion fk{Ek) = OkEk and the proportional-fair- 
energy (PFE) maximization where fk{Ek) = In(ii^fc). 

V. Optimal Solutions for Resource Allocation 

In this seetion, we obtain the optimal solutions for resouree alloeation for WET in a baekseatter 
eommunieation system, for the single-ER ease and the multiple-ER ease in Seetion IV-AI and Seetion 
IV-Bl respeetively. 


A. Optimal Solution for WET to Single ER 

In this seetion, we analyze the optimal solution for WET to single ER, i.e., K 
From (fTSl) . the harvested energy by the ER is given by 


Efq) = M(3i{p^yfr - q) 


1 - 


ct 2(M- 1) In (1 + 


Ph 


Piq 


= i,ei = 1- 


( 22 ) 


We note that it suffiees to maximize Ei{q), subjeet to 0 < g < PayfT, sinee the utility funetion 
fi{Ei{q)) is assumed to be a monotonieally inereasing funetion of Ei{q), for both WSE maximiza¬ 
tion and PFE maximization. We obtain the optimal solution to Problem (Pi) for WET to single ER, 
in the following theorem. 


Theorem 1 (Optimal Solution for WET to Single ER). The optimal solution to Problem (Pi) for 
WET to single ER is given by 


/3irpave(M—1) 


q^ = 


0 , > M 

g*, otherwise 


(23) 


where q\ is the unique solution for q G (0,Pave7"] to the equation 
fM-l) f /3fq , f /3fq , \ /)iV /SiV 


Pa 


M 


lMa2 ^ 




MTa2 


/3f{M - l)paveg 
MV2 


= 0 . 


(24) 


Proof: (Sketeh) When the noise varianee the objeetive funetion Ei{q) ean 

be easily shown to be monotonieally deereasing with respeet to g. The optimal solution is thus zero. 
When cr^ < > the objeetive funetion Ei{q) ean be shown to be strietly eoneave with 

respeet to g for g G (0,pave7']. The solution is thus unique. See details in Appendix O ■ 
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Theorem [T] implies that when the noise level for CE at the ET is too high, it is better for the ET 
to broadeast energy in all direetions, without beamforming. This is as expeeted, sinee the energy 
used for CE does not justify the gain aehieved from beamforming. 


B. Optimal Solution for WET to Multiple ERs 


In this seetion, we maximize the total utility for WET to multiple ERs. In partieular, Seetion IV-B 1 1 
and Seetion IV-B2I eonsider the WSE maximization and the PEE maximization, respeetively. 

1) Weighted-Sum-Energy Maximization: Eor WSE maximization, the total utility is rewritten as 

K K 


u{q, 0 = {p^yeT -q)J2 Mk + ak{q)^k, 

k=l k=l 


(25) 


where the funetion ak{q) is given by 


ak{q) = 0k/3k{M - l)(paveT - q) 



KMa^ 


In 



Y 


(26) 


Erom the inequality ln(a; + 1) < x, Vx > 0, it is obvious that ak{q) > 0 for g > 0. It is noted that 
afc(g) is a strietly concave function of g G f),PaMeT). 

We obtain the optimal solution to the WSE maximization problem in the following theorem. 


Theorem 2. The optimal energy allocation weights are 


= 


1, k = k* 

0, otherwise 


(27) 


where k* is given by the following criterion 

k* = argmax ak,{qlj + {paveT - qljYl 


K 


ki 


k=l 


where q^^ is the unique solution that maximizes the function ak^{q) for q G (0,Pave7"]- The corre¬ 
sponding optimal training energy g* = g^*. 


Proof: Erom (1251) . for any given g G (0,pave7'], the objective function is a linear function of 
.^fc’s. To maximize the weighted sum of harvested energy, it suffices to allocate all energy to single 
energy beam toward the ER with the largest linear-combination weight afc(g). This gives (1271) . Eor 
g G (0,pave7'], the solution that maximizes the function afc^(g) is unique, as ak{q) is strictly concave. 
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This completes the proof. ■ 

Theorem [2] implies that only one beam is used to transfer energy to one partieular ER in the WET 
phase, although the ET has ehannel estimates for all ERs. Energy is thus wasted for estimating 
ehannels of other ERs. Moreover, when more antennas are employed at the ET, the only energy 
beam beeomes more eoneentrated. Other ERs without dedieated energy beam ean harvested very 
little energy, resulting in severe unfairness among ERs. This observation will be numerieally verified 
in Seetion IVIIi Henee, we will eonsider another utility funetion that takes fairness among ERs into 
eonsideration, in the next seetion. 

2) Proportional-Fair-Energy Maximization: In this seetion, we aim to maximize the log-sum of 
the energy harvested by all ERs, whieh is known to result in proportional fairness ifTTl . Erom (flSl) . 
the total utility is rewritten as 

K 

U{q,$.) = '^ln{bk{q)^k + dk{q)) , (28) 

k=l 


where the quantities depending on q are given by 


h{.q) = /3k{M - l){p^^eT - q) 



KMa^ 


In 



KMa^ )\ ’ 


4(g) = 4 {PaveT - q ). 


(29) 

(30) 


Henee, the PEE maximization problem is rewritten as 

K 


(P 2 ) min - In (4(g)4 + 4(g)) 

(31a) 

s. t. eonstraint (|21b|). (I21el). (|21d|) 

(31b) 


Before further analysis, we give the following definitions IfTSl . 

Definition 1. A function g : X xy ^ TZ is called biconvex, if g{x, y) is convex in y for fixed x E X 
and is convex in x for fixed y E y. 

Definition 2. A problem is a biconvex problem, if it optimizes a biconvex function over a given 
biconvex or compact set. 

We then have the following theorem for the Problem (P 2 ). 
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Theorem 3. The Problem (P 2 ) is a biconvex problem. 


Proof: We note that the logarithm funetion is eoneave and inereasing. Given q, the summation 
term in (I31al) is eoneave, as it is a eomposition of a eoneave and inereasing funetion (i.e., ln(-)) and 
a eoneave funetion (i.e., linear funetion of ^k)- The objeetive funetion (I31al) is thus strietly eonvex. 

On the other hand, given the argument of the logarithm funetion is eoneave, as bk{q) is eoneave 
(see Appendix O- The objeetive funetion (I31al) is also strietly eonvex, as it is the sum of a family 
of eompositions of a eonvex and deereasing funetion (i.e., — ln(-)) and a eoneave funetion (i.e., 
linear funetion of Ek{q)). Clearly, the domain for Problem (P 2 ) is a eonvex set. By Definition [2l 
the Problem (P 2 ) is a bieonvex problem. ■ 

In general, a bieonvex problem is noneonvex and has multiple optima. Before giving the algorithm 
to find solution for the bieonvex Problem (P 2 ), we first deeompose the problem into two subproblems. 
For Subproblem (P 2 a), given the training energy q, we optimize the energy alloeation weights 
namely, this performs beamforming energy alloeation. For Subproblem (P 2 b), given the energy 
alloeation weights we optimize the training energy q, namely, this performs training energy 
alloeation. 

Subproblem (P 2 a) is a (strietly) eonvex optimization problem. There is thus a unique global 
solution. We give the optimal beamforming energy alloeation for Subproblem (P 2 a )7 in Theorem |4l 


Theorem 4 (Optimal Beamforming Energy Alloeation). Given training energy q G (0,paveT], the 
optimal solution for subproblem (P 2 a) follows a water-filling method. In particular, the optimal 
beamforming energy allocation weight is given by 


Ckiq) 


max 



4(g) 1 
4(g) J’ 


where the water-level is the unique solution to the equation 



dkjq) \ 
bk{q) J 


1 . 


(32) 


(33) 


Proof: It is proved by using Karush-Kuhn-Tueker (KKT) eonditions. See Appendix iDl ■ 
As shown in the proof for Theorem [3l Subproblem (P 2 b) strietly eonvex. We ean thus find the 
unique optimal solution q\ e (0,PaveT], by using any eonvex optimization toolbox, sueh as [|T9ll . 
although it is diffieult to obtain the elosed-form solution. We note that the objeetive funetion (I31al) 
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Algorithm 1 Algorithm for PFE Maximization: 

1 : Parameters: {/3fc}, Pave, T, M, K, cr^ e. 

2 : Initialization: Choose any go ^ {0,PaveT]), and any feasible that satisfies (I21cl) and (I21dl) . 

compute Uo = t/(go, set Ui = Uo + 2e, t = 0. 

3: while \Ut+i — Ut\ > e do 

4: Keep qt fixed, use the water-filling results in Theorem |4] to find the optimal energy allocation 

weights ^t+i- 

5: Keep ^t+i fixed, find qt+i that minimizes the objective function in (I31al) . by using standard 

convex optimization techniques. 

6: t = t + l. 

7: Use (l28l) to compute the updated utility as Ut+i = U{qt, $,t)- 

8 : end while 

9: return g* = qt, t^max = U{qt,^t)- 


is derived for WET via energy beamforming, and thus not applicable for the case of g = 0 in which 
omnidirectional transmission is used. Eor that case, the total utility is a constant (PkPaveT), 

regardless of the choose of .^^’s. Eor the case of g = Pa^eT, the total utility is — cxd, as no energy 
is harvested (i.e., no time is allocated for WET). Hence, the final optimal g* for given ^ is chosen 
between g* and zero. 

Eor a biconvex problem, there is no algorithm that ensures to find the global optima [|T^ . In the 
state-of-the-art literature, the block coordinate descent (BCD) algorithm is computationally efficient 
and with performance guarantee, as it ensures to converge to a partial optimal solution [flOll . Hence, 
we propose a BCD-based Algorithm [T] as follows. 

VI. Discussion 

In this section, we discuss some extension work and interesting practical issues. Eirst, the analysis 
and results in this paper are also applicable to WET in a time-division-multiplexing (TDM) manner. 
The WET phase of (T — r) symbol periods is divided into K slots. The k-\h WET slot consists 
of ^fc(T — r) successive symbol periods. The relative time allocation coefficients are subject 
to Ef= = 1- In the fc-th WET slot, the ET delivers energy to ER k via energy beamforming 
using only the estimated BS-CSI ak for ER k. It can be shown that for WET in a TDM manner, 
the harvested energy is the same as in (fT71) . and all the resource allocation results are also the same 
as those for the current energy multicasting in this paper. Hence, the WET via energy multicasting 
is equivalent to WET in a TDM manner. This is basically because the amount of harvested energy 
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has a linear relationship with either the energy-harvesting time or the power of the energy beam. In 
this paper, we use multiple beams to transfer energy to all ERs eoneurrently, as this operation ean 
avoid frequent switehes between WET to different ERs, in eaeh frame. 

Based on current work, other issues can be further addressed, such as WET under nonlinear 
backscattering and channel reciprocity. Eor ERs with nonlinear backscattering, the reflection coef¬ 
ficient pk varies with the incident power as a backscatter ER is typically matched at some 
fixed power level, and thus mismatched at other incident power levels. Given the ER hardware, the 
training energy can also be optimized, once the reflection characteristic is experimentally measured. 

On the other hand, to achieve smaller size and lower cost, some REID readers use single antenna 
for both transmission and reception, by introducing an RE isolator such as a circulator or a directional 
coupler [l22l . With the assumption of reciprocal channels, in ll23l . the phase of the forward channel 
of each reader antenna was estimated. However, there may be phase ambiguity of integer multiple 
of TT, which may result into less power delivered to the RE tag. Nonlinear methods may achieve 
better estimation of forward channels, but are beyond the scope of this paper. 

VII. Numerical Results 

In this section, numerical simulations are given to corroborate our analysis. We assume that the 
average power is Pave = 1 W. We set the number of transmit antennas as M = 4. We assume each 
frame consists of T = 200 symbol periods. Eor convenience, we normalize each symbol period to 
be one second, resulting into 200 J of energy consumption in each frame. The carrier frequency is 
5 GHz, and the bandwidth is 100 KHz. We set the power spectrum density of noise as —140 dBm/Hz, 
which implies the noise power at the ET is = —90 dBm. We consider two ERs, i.e., K = 2. 
We take the path loss model as 10“^/!“^, where the path loss exponent is 3, and D is the distance 
between the ET and an ER. A 30dB path loss is assumed at a reference distance of 1 m. The energy 
harvesting efficiency at each ER is assumed to be p = 0.8. The reflection coefficient for both ERs 
is Pi = p 2 = 0.8 -f 0.5i. All the simulation is based on 100, 000 Monte Carlo simulation runs. 

A. Single-ER Case 

We first consider the single-ER case. We assume the distance Di = 6 m, which implies the path 
loss is /9i = 4.6296 x 10“®. We first simulate the harvested energy for two benchmarks, i.e., the 
case of the perfect E-CSI and the case of no CSI. Eor both cases, no channel-training time nor 
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energy is required. With the perfeet F-CSI, the ET performs maximum-ratio-transmit (MRT), and 
the harvested energy is obtained as 2.96 mJ. For the ease of no CSI, the ET performs omnidireetional 
transmission without beamforming, and the harvested energy is 0.74 mJ. 

We then simulate the harvested energy by using the estimated BS-CSF Fig. [3] plots the harvested 
energy versus the training energy q. We observe that when the training energy approaehes zero, 
the harvested energy is 0.75 mJ, whieh approaehes that for omnidireetional transmission. This is 
beeause no BS-CSI ean be inferred, due to zero training energy. 

By simulation, the optimal training energy is ^ = 9.0 J. The maximally harvested energy is 2.68 
mJ. From Theorem [H the optimal training energy is g* = 9.2 J, whieh eorroborates the simulation 
results. More importantly, we observe that eompared to the MRT seheme, the maximally harvested 
energy via energy beamforming by using the estimated BS-CSI suffers only a slight reduetion 
of 9.6%. As expeeted, the harvested energy is also inereased signifieantly by 262%, eompared 
to omnidireetional transmission. On the other hand, we observe that the dash-dot o-marker eurve 
obtained by analysis eoineides with the solid o-marker eurve obtained by simulations. We also see 
that the dashed D-marker eurve obtained by the lower bound is tight, whieh is obtained analytieally. 
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B. Comparison of Harvested Energy 

In this section, we compare the harvested energy to the net harvested energy by using the estimated 
F-CSI in a traditional radio communication system, to the harvested energy by using omnidirectional 
transmission, and to the harvested energy by using the perfect F-CSI. We assume that the traditional 
radio communication system operates in time-division-duplex mode. Thus, the FT has to first send 
forward pilots, such that the ER can estimate and feed back the F-CSI; additionally, the ET has to 
estimate backward channels by receiving pilots sent from the ER, and then recover the estimated 
E-CSI. Eor the ER, we consider the energy used for sending backward pilots and feeding back the 
estimated E-CSI, neglecting the energy used for data acquisition and computation. In particular, the 
ER employs analog feedback [12^ . as it requires a lower feedback rate and has a smaller feedback 
delay. Following the scheme in Section II of lt2^ . the ET performs minimum-mean-square-error 
(MMSE) estimation for the backward channels, and uses the optimal MMSE filter to recover the 
estimated F-CSI. We denote the power for backward transmissions by p^. For the ER, the total 
energy used for backward transmissions is Eq = {M + l)pu- In order to maximize the net harvested 
energy (after subtracting Eq), there is an optimal pu- In the sequel, the maximally net harvested 
energy is obtained by jointly optimizing both the forward resource allocation and p^. 

Fig. m compares the harvested energy and the efficiency of different schemes in the upper and 
the lower figure, respectively. The efficiency is defined as the ratio of the harvested energy divided 
by the harvested energy by using the perfect F-CSI. As expected, the efficiency for omnidirectional 
transmission is 0.25, as the harvested energy via energy beamforming by using the perfect F-CSI 
is increased by M = 4 times, compared to that by using omnidirectional transmissions. 

We observe that when the average transmit power pave at the ET exceeds 4 W, energy beamforming 
by using the estimated F-CSI and the estimated BS-CSI achieve the same efficiencjy of about 96%. 
For smaller pave, the efficiency of energy beamforming by using the estimated BS-CSI still exceeds 
90%, and is slightly lower than that by using the estimated F-CSI. For instance, for pave = 2 W, the 
harvested energy by using the estimated BS-CSI (i.e., 5.5 ml) suffers slight reduction, compared 
to that (i.e., 5.7 ml) by using the estimated F-CSI. After normalizing to the harvested energy (i.e., 
5.92 ml) by using the perfect F-CSI, the efficiency of energy beamforming by using the estimated 


''We observe that the efficiency by using the estimated F-CSI is almost constant. This is because for each pave at the ET, the 
estimated F-CSI at the ER is near perfect due to only one-way forward-channel propagation, and the optimal energy used for 
backward transmissions is numerically shown to be the same. 
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BS-CSI is 93%, which is slightly degraded, eompared to the effieieney (i.e., 96%) by using the 
estimated F-CSI. This is eneouraging, as almost all eomplexity of hardware and eomputation at 
the ERs is shifted to the ET, at the eost of slight reduetion in the harvested energy. The energy 
beamforming by using the estimated BS-CSI is thus effieient and attraetive for transferring energy 
to (ultra-)low-power and low-eost wireless deviees. 



Average ET transmit power (W) 


^ 0.8 
•G 0.6 
H 0.4 
0.2 


^ gggooooogooo -i i 




E.B. using BS-CSI 
^— E.B. using F-CSI 
— Omni, trans. without E.B 


AAAAAAAAAAAAAAAAAAA. i 


2 4 6 8 

Average ET transmit power (W) 


10 


Fig. 4: Harvested energy and effieieney for different energy beamforming (E.B.) sehemes 


C. Multiple-ER Case 

In this seetion, we simulate both the WSE and the PFE maximization problems, for the two-ER 
ease. Eor WSE maximization, we fix Di = 4 m and ^2 = 6 m, and ehoose the eombination weights 
9i = 0.3 and 62 = 0.7, to balanee the energy harvested by the nearer ER 1 and the further ER 2. 
Erom simulation, the optimal training energy is g* = 4.1 J and the optimal weight is = 1. The 
harvested energy is 9.71 ml and 0.7431 ml for ER 1 and for ER 2, respeetively. Erom Theorem [2l 
we have q* = 4.0, = 1. The harvested energy is 9.66 ml and 0.739 ml for ER 1 and for ER 2, 

respeetively. The simulation results eorroborate the analysis. 

Eor PFE maximization, we fix Hi = 4 m, and ehoose different distanees D 2 > Di. The results 
are given in Table HI in whieh we use the normal notation (e.g., q and ^k) for the analytie results. 
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and use the notations with hat (g and respectively) for the numerical results. We observe that the 
obtained solution by using the BCD-based algorithm is close to the optimal solution. The simulations 
corroborate the analysis. 


TABLE I: Optimal solutions versus distances (PEE Maximization) 


D2 

q* 

g*(BCD) 

r 


?f(BCD) 


El 

El 

E 2 

E 2 

4 

3.83 

3.80 

3.83 

0.5000 

0.4998 

0.5001 

6.052E-3 

6.076E-3 

5.908E-3 

5.932E-3 

5 

5.46 

5.42 

5.45 

0.5112 

0.5105 

0.5102 

6.040E-3 

6.071E-3 

2.915E-3 

2.93 lE-3 

6 

7.45 

7.43 

7.44 

0.5285 

0.5265 

0.5248 

6.148E-3 

6.169E-3 

1.563E-3 

1.581E-3 

7 

9.65 

9.63 

9.64 

0.5502 

0.5485 

0.5476 

6.238E-3 

6.249E-3 

9.006E-4 

9.019E-4 

8 

10.86 

10.84 

10.85 

0.5898 

0.5876 

0.5862 

6.490E-3 

6.507E-3 

5.299E-4 

5.302E-4 


Moreover, we compare the harvested energy for the WSE maximization problem and the PEE 
maximization problem, for different number of antennas M deployed at the ET. We fix Di = 
4 m, ^2 = 6 m. Eig. [5] plots the maximally harvested energy of each ER versus M. Eor WSE 
maximization, we observe that as M increases, the harvested energy by the nearer ER 1 increases, 
while the harvested energy by the further ER 2 remains as a small constant. This is because the 
ET uses only one energy beam toward the ER 1. Therefore, the harvested energy is unfair among 
ERs. Eor PEE maximization, however, the harvested energy by both ERs increases, as M increases. 
Compared to WSE maximization, the harvested energy by the further ER 2 is increased significantly, 
although the energy harvested by the near ER 1 is about half of that for WSE maximization. The 
harvested energy is more balanced between the two ERs. Hence, we conclude that better fairness 
is achieved by PEE maximization. 


VIII. Conclusion 

This paper studies the energy beamforming by using the estimated BS-CSI. We analyze the 
harvested energy, by investigating the effect of the ambiguity of backward channels. Moreover, we 
optimize the channel-training energy and the energy allocation weights, for two utility maximization 
problems. Eor WSE maximization, the optimal WET scheme is to use one energy beam, resulting 
in unfairness of harvested energy. Eor PEE maximization, we show that the problem is biconvex, 
and propose a BCD-based algorithm to find the close-to-optimal solution. The harvested energy 
by using the estimated BS-CSI is numerically shown to suffer slight reduction, compared to that 
by using the perfect E-CSI, and also to the net harvested energy by using the estimated E-CSI in 
traditional radio communication systems. Hence, the energy beamforming by using the estimated 
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Fig. 5: Maximally harvested energy versus M. 


BS-CSI is a promising WET seheme, espeeially for transferring energy to (ultra-)low-power and 
low-cost wireless devices that are neither capable of estimating the channels nor sending pilots 
actively. Other interesting issues are also discussed and remain to be addressed, such as the WET 
under nonlinear backscattering and channel reciprocity. 

Appendix A 
Proof for Eemma[I] 

Recall the harvested energy given in ([7]). The term in the round brackets of ([7]) is the sum of a 
scaled identity matrix and a rank-one matrix. The eigenvectors can be constructed as follows: take 
the normalized as the right eigenvector corresponding to the maximal eigenvalue, and construct 
other mutually orthogonal eigenvectors by Gram-Schmidt algorithm. The term in the round brackets 
of ([7]) is then maximized, when the beamformer is given by 

Wfc(afc) = (34) 

||hfc|| 

Clearly, the argument of the inner expectation is still maximized, when a common angle is 
introduced to all beamforming weights. Hence, an optimal beamformer that depends on the available 
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estimate of the baekseatter ehannel is given by 


Wfc(afc) = 


Jhfcll ll^fcll 

With the beamformer in (1351) . the maximally harvested energy in ([7]) is rewritten by 


( 35 ) 


Eo,k(,r,pi) = P2E, 


9k 


Pk(jlk{9k) ^ 


hfh, 


[Pk + crlM + 

1 


= 


f^k + (^lki9k) 


= p2Mf3k f 1 - 


'9k 


/3fcTpi|gfcP I 1 
KMa'^ '■ 


(36) 


where (a) is obtained by substituting E 


hfelgfe 


= M (/Sfc + crl^kipk)), whieh is from the faet 
that the estimated F-CSI is eonditionally distributed as hk\gk ~ CJ\f (Om, {Pk + cre,ki9k)) ^m), (b) 
is from the eonditional error varianee given by dH). 


Appendix B 
Proof of Lemma [2] 

Substituting (fTTI) into (fT2l) . the harvested energy by ER k is rewritten as 


Ek{T,pi,p2,^k) = 
P2{T-t) 



( 

ik hf afc f 

't) 


l|afc|P 


V 



+ E^Ej^^|_^ 


1/^ (hf^afc) ^ y/^hjkat 


afc 2 


i^k 


a* 2 


(37) 




\f^^k ^k (h^a,) 


H 


afc 2 


i^k 


a* 2 


+ ®^A®^hJA 




EE 

jj^k ij^k 


a, 2 ah 2 


\ 

J/ 


In the sequel, we investigate the four terms in the outer round braeket of (l37l) . Given g^, the 
random variable + e^. Reeall that eonditioned on a^ and gk, the distribution of is 

9k 

given by dH). Conditioned on the ehannel estimate A, the first term is rewritten as 


EjE. \2 

A hfc A 


ik |h^afc| 

llafclP 




a^ ^hk\ak,9k 


[hfchf ] fifc 


^k 


= MMk 



M -l 
M 


E 


'9k 


PfcT-pilgfcp I 1 
. AM(t 2 f -L 


(38) 
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where (a) is obtained by following similar steps in the proof (Appendix for Lemma [B 
The seeond term in (lT7l) is rewritten as 
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EtE, IX 
A hfe A 




afc 2 


i^k 


^* 2 


i^k 


ll^fcll 11^*11 


(a) 


9 fe ,ai I gt, 

i^k 


—;—9“7^—r-'-AT + :—^—;-7 


afcll \Pk + crlkigk) |5(fc|2 + al^{gk)y 




(39) 


where (a) is from ([3]), and (b) is duo to the faet that and ||5A are independent zero-mean 

ll&fcll ll^fll 

random vectors, for any i ^ k. 

The third term in (ITTl) . which is the conjugate of the second term in (1371) . is similarly obtained as 


®"A®hfc|A 




H 


afc 2 


i^k 


ad 2 


= 0 . 


(40) 


The fourth term in (l37l) is rewritten as 




VMi (hf%)^hfa, 


EE 

j¥^k l^k 


(a) 


(fe) 


(0 


af ( E,, 


0Eas.,aj 

'CjEsj.,sj 

j¥=k 

= /3fc(l - ^fc) 

3¥^k 


ad 2 a; 2 


!? (Eh,|S. [h,h«]) 


EE 

j^k l^k 


if (Eh,|Sfe [hfchf])a; 

11 a? II 11 a/11 


/3feO-e fe(9fc) 


Im + 




yfc+<^,,A9fc) |9fcP(/3fc+<fc(9fc)) 


(41) 


where (a) is from the fact that -py and are independent zero-mean random vectors, for any 
j 7 ^ /, (b) is from ([3]), (c) is from the fact that conditioned on gk, the vector is distributed as 
CM (Om, Is'fcP {f3k + o'e,fc(5'fc)))’ ^ad (d) comes from the normalization condition = 3- 

Substituting (1381) . (I39l) . (I40l) and (l4T]) into (ITT]) , we obtain the harvested energy as in (fT3l) . 
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Appendix C 
Proof of Theorem [H 

For convenience, the objective function in (l22l) is rewritten as 

C 2 ln(a; + 1) 


Ei{x) = M/3iT(pave - CiX) 


1 - 


X 


( 42 ) 


where Ci = the introduced variable x = G (0,PaveCi and C 2 = G (0,1). 
The first derivative is derived as follows 


E[{x) = 


Mf3iT 


[c2Pave(a; + 1) ln(a; + 1) - CiX^ - Ci(l - C2)a;^ - C2Pavea^] • 


x‘^{x + 1 ) 

From (1431) . the second derivative is then derived as follows 


(43) 


E'lix) = 


C2PaveMl3iT 

x‘^{x + 1 ) 

M C2PaveM/)iT 

x‘^{x + 1 ) 

= -C 2 PaveM/)iT 


2(x + 1) ln(a: + 1) x 


X 


a: + 1 




Pa 


2(a: + l) 1 

'J' 'Yl 


E 

n=3 


X 


n=l 
n—3 


X 


n Vx + 1 

+ 


X 


n{x + 1)^ Pave(a; + 1)^ 


X + 1 
< 0 


CiX 

Pave 


(44) 


where (a) is from the Taylor’s expansion ln(x + 1) = ^ ^ “|- 

From L'Hopital’s rule, the right limit at x = 0 for the first derivative function is given by 

hm E[(x) = MI3X lim + P = [c.p.„- 2 ci(l - c.)] 

x^o+ ^ ^ x^o+ 6 x + 2 2 ^ ^ 


We consider two scenarios. First, when < 


/3lXpa„e(Vf-l) 

M 


(i.e.. Cl < 2 ( 1 ?^)’ *S 2 ]), we have 

lim E[{x) > 0. There thus exists a unique solution x* G (0,paveCj“^], since i?i(x) is concave with 

x—>- 0 + ^ 

respect to x and ii^j(pavecr^) = 0 - 


Second, when > 


/ 3 irpave(M — 1 ) 


(i.e.. Cl > 5773 ^), from (1451) . we have lim E[{x) < 0. 


M ^ 2 { 1 -C 2 )^ — ^^ 0 + 

Moreover, from the Taylor’s expansion of ln(x + 1), the first derivative in (l43l) is rewritten as 


E'Ax) = 


Mf3iT 
X + 1 


C2Pa 


W M/)iTx 
“ X + 1 


2(x + l) 

C 2 Pave ^ ^ 


+ C 2 Pave “ 


°o ^ ^n-2 


n (x + 1)’ 

n=3 


— - CiX - Ci(l - C2) 


°o 1 ^n-3 


^ n (x + l)”“i 


- Cl 


n=3 
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W M/3iTci , , 
< —- ^9{x) 


(46) 


x{x + 1)' 

where (a) is from the faet C 2 Pave < ci(l — C 2 ) and a; > 0, (b) is due to the faet C 2 Pave < ci(l — C 2 ), 
and the function g{x) is given by 


9{x) = 


2 ( 1 -C 2 )X) 


^n—1 


X 


^n{x + 


X 


n=3 


(47) 


Clearly, we have g{0) = 0. The derivative of g{x) is derived as 


9'{x) = 


2(1 — C2)x (n — l)a; 




n—1 


x+l ^n{x + lY~^ 

n=3 ^ ' 


2x 


(a) 4(1 _ C 2 ) “ 


X 


n—1 


- (X + 


— 2x = X 


4(1 - C 2 ) 
3(a; + 1)2 


- 2 


(b) 

< 0 


(48) 


where (a) is from x > 0, and (b) is from 0 < C 2 < 1 and x > 0. The derivative in (1461) is always 
non-positive, due to g{x) < 0, Vx G (0,PaveC)"^] ■ Hence, when cr2 > ^ optimal 

X* = 0. For ^ optimal x* is the unique solution for E[{x) = 0, where the first 

derivative is given by (|4^ . This completes the proof. 


Appendix D 
Proof for Theorem [4] 

For notational convenience, we omit the given argument q in this proof. The Lagrangian is 
constructed as follows 

K K / K 


• • • , Ai, • • • , Xk, b') = — ^ \n{bk^k + dk) — ^ Xk^k + b' j ^ ■Cfc — 1 


k=l k=l 

The Karush-Kuhn-Tucker (KKT) conditions are thus written as 


,fc=i 


(49) 


Ck > 0, Vfc 


K 




k=l 

i * 


> 0, Wk 


KCk = 0 


bk^t + dk 


- A^ + I/* = 0, Vfc 


(50) 

(51) 

(52) 

(53) 

(54) 
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where the equality (a) is derived by taking the first derivative of the Lagrangian. We note that 
is a slack variable which can be eliminated. The conditions (l52l) (15^ (l54l) are thus rewritten as 


u* > 


u 


+ dk 

h 


= 0 


bk^l + dk 


Wk 


(55) 

(56) 


If I/* < the condition (l56l) can only hold for > 0. Then the condition (l55l) implies 
= bkl*+dk • Equivalently, we have ■Cfc = ^ ~ E then > 0 is impossible. This is 

because the complementary slackness condition (l55l) is violated, as z/ > , if Q > 0. That is, 

we have = 0, when u* > ^. Hence, we obtain the optimal WET time allocation coefficients 


Ck = max <! 0, 4 ^ • 


z/’ 


Substituting (l57l) into the condition (ISTl) . we obtain 

K 


Vmaxjo.T-^l 


= 1 . 


(57) 


(58) 


The left hand side of (1581) is a piece-wise linear increasing function of with breaking point at 

The equation (15^ thus has a unique solution that is readily determined. This completes the proof. 
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